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Abstract: 

We compute the t — ► WZ b decay rate, in the Standard Model, at the leading order in 
perturbation theory, with special attention to the effects of the finite widths of the W and 
Z bosons. These effects are extremely important, since the t —>■ WZb decay occurs near 
its kinematical threshold. They increase the value of the decay rate by orders of magnitude 
near threshold or allow it below the nominal threshold. We discuss a procedure to take 
into account the finite-width effects and compare the results with previous studies of this 
decay. Within the Standard Model, for a top quark mass in the range between 170 and 
180 GeV, we find BR(t -> WZb) ~ 2 x 10" 6 , which makes the observation at the LHC 
very difficult if at all possible. 



1 Introduction 



The t — > WZb process is a particularly interesting rare decay mode. In the Standard 
Model, all rare top decays are extremely suppressed with respect to the largely dominant 
t — ► Wb decay. As possibly it could be observable at the LHC, the t — > WZb decay has 
already attracted considerable attention in the literature 0-[||]. 

The most peculiar feature of the t — ► WZb decay is the fact that the process occurs 
near the kinematical threshold, m t ~ my/ + rriz + As observed in ref. 0, it is then 
crucial, in the theoretical evaluation of the decay rate, to take into account the finite-width 
effects for the W and Z bosons. For a top quark mass around approximately 175 GeV, 
these effects increase the value of the decay rate by some orders of magnitude. 

Even with the sizeable enhancement of the rate induced by finite-width effects, the 
branching ratio for t —>■ WZb is predicted to be of the order of 10~ 6 , which could be too 
small for the decay to be observed even at the LHC, where 10 7 — 10 8 top quark pairs are 
expected to be produced per year. Therefore, the observation of this decay mode at the 
LHC could be a signal of physics beyond the Standard Model (see, for example, ref. [|]. ||). 

In this paper we first discuss the results for the t —>■ WZ b decay rate computed at tree- 
level, in the Standard Model, by neglecting the effects of the finite widths of the W and Z 
bosons, i.e. by treating these particles as stable particles in the final state. This quantity 
has also been computed in refs. [0J-|| but the results of these studies are in disagreement. 
Our results, in the limit of vanishing widths, will be presented in sec. 0. They agree with 
those obtained in ref. Q, but present numerical differences with respect to those of ref. 
and, near threshold, of ref. [0. 

In sec. H], we compute the t —>■ WZb decay rate by taking into account, as in ref. ||, 
the finite-width effects of the W and Z bosons. We consider two different approaches. The 
first approach, which we will refer to as the convolution method, has been discussed, in a 
different context, in ref. M. It is based on a convolution of the t — ► WZb decay rate with 
two Breit-Wigner-like distributions for the invariant masses of the W and Z bosons. These 
distributions correspond to the imaginary part of the gauge boson propagators, and are 
centered around the physical values of the W and Z masses with spread controlled by the 
physical gauge boson widths. 

The second approach has been followed in ref. [Qj. It consists in considering the decay 
rate for the process t — > W*Z*b followed by the decays of the virtual W and Z bosons, e.g. 
W* — > fiv^ and Z* — > e + e~ . The t — > WZ b decay rate is then obtained by dividing the rate 
computed for the full decay chain by the product of the W — > [w^ and Z — > e + e~ branching 
ratios. We will refer to this second approach as the decay-chain method. Incidentally, 
we note that one Feynman diagram entering at tree-level has not been included in the 
calculation of ref. f2|. 

A difficulty inherent to the decay-chain method is that other initial processes, besides 
the t —>■ WZb decay of interest, may contribute to the decay chain. For instance, when 
an electron pair is considered in the final state, a virtual photon, instead of the Z boson, 
may be produced in the intermediate states. Therefore, in this case, the direct connection 
between the full t — > bybV^e' process (total number of events) and the t — > WZ b decay 
("signal") is lost. In order to suppress the contribution of the virtual photon, a kinematical 



1 



z 



b 



W 



b 



b 



W 




< 




< 




t 



t 



t 



b 



t 



W 



Z 



z 



Figure 1: Feynman diagrams contributing at tree-level to the t — > WZb decay. 



cut on the invariant mass of the electron pair (requiring m ee > 0.8 mz) has been introduced 
in ref. f2|. This cut is indeed effective and, in the selected region of the phase space, the 
t — > h[iv tl e + e~ process mostly proceeds through the initial t —>■ WZ b decay. However, with 
this requirement, the definition itself of the t — > WZb decay rate becomes dependent on 
the specific choice of the cut. The purpose of this paper is to show that, by using the 
convolution method, a convenient definition of the t — > WZ b decay rate can be provided, 
which is independent of any kinematical cut. In this way, we find results for T(t — > WZb) 
which differ by approximately a factor 3 with respect to those obtained following the 
procedure of ref. [[J. Of course, the full calculation of t — > fe/iz/^e+e - , and the analysis 
of the kinematical cuts, can be relevant for the evaluation of the measurable signal plus 
background rate. 

The calculation of the QCD radiative corrections goes beyond the scope of the present 
paper. However their effect could be important and give rise to a further reduction of the 
rate due to Sudakov suppression factors that tend to become dominant near the end of 
phase space. 

2 The t — > WZb decay rate in the limit of stable W 
and Z bosons 

Within the Standard Model, three Feynman diagrams contribute to the t — > WZ b decay at 
the leading order in perturbation theory. These diagrams are shown in Fig. |l] and represent 
the amplitudes with the final Z boson radiated either by the initial top, or by the final b 
or W respectively. 

We have computed the decay rate for this process by considering first the final W and 
Z bosons as stable particles, i.e. neglecting their finite widths. In table [I] we collect the 
values of the various quantities which have been used to obtain all the results presented 
in this paper. The integration over the final phase space of the three massive particles 
has been performed numerically. For the values of masses considered in this paper, the 
kinematical threshold for the process is at m t = 176 GeV. Note that the exact threshold 
value depends on m^. At the parton level, rrib is in principle scale dependent. A low scale 
appears justified here because of both the very limited phase space allowed to a virtual b 



2 



m b 


mw 


m z 


sin 2 6V 




4.5 GeV 


80.3 GeV 


91.2 GeV 


0.23 


1.166 x 10~ 5 GeV" 2 


IV 


Tz 


BR(W -> /i^) 


-> e + e~) 


BR(Z -> i/p) 


2.06 GeV 


2.49 GeV 


0.102 


0.03367 


0.06667 



Table 1: Values of the numerical constants used in this paper 

and the value of lightest 5-meson mass. 

In the limit of neglecting the finite widths of the W and Z bosons, the decay is forbidden 
below the threshold. For larger values of the top quark mass, we obtain the values of the 
branching ratio presented in table [|[] The uncertainties on these results, coming from the 
numerical integrations, are estimated to be smaller than the last digit shown in the table. 
For values of the top quark mass around 180 GeV, we find that the branching ratio is 
extremely small, of the order of 10~ 8 . 

The results obtained in this section agree with those of ref. ||. We also agree with the 
analytical expression of the tree-level Feynman amplitude squared published in that paper. 
On the other hand, we disagree with the numerical results of ref. PJ (we note that there 
is an apparent error in the integration limits over the phase space of the three massive 
particles (eq.(40) of that paper)). 

We also found discrepancies, of the order of 50%, with respect to the results obtained in 
ref. [|[ in what they call the narrow-width approximation. This procedure consists in deri- 
ving the values of the t — > WZ b decay rate starting from the calculation of V(t — > bfiue + e~) 
and constraining the invariant masses of the virtual W and Z bosons to their physical 
values. In turn, these constraints are achieved by introducing ad hoc delta functions into 
the phase space integrals. We find that, in the stable particles limit, our results differ from 
those of ref. by a factor 2 for m top = 177 GeV and by approximately 10% for m top = 185 
GeV. 



throughout this paper, we approximate the total decay width of the top quark with the partial rate 
r(f — > Wb), computed at the corresponding values of the top quark mass. Therefore, we define for instance 
the t — » WZ b branching ratio as: 
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BR(t -> WZb) 
Stable Particles 


BR(t -> WZb) 
Conv. Meth. 


BR(t — > bfiu^uu) 
BR(W -> [iv v )BR{Z -> vv) 


150 


- 


0.65(3) • 10~ 6 


0.79(1) • 10~ 6 


155 


- 


0.90(5) • 10~ 6 


0.87(1) • 10~ 6 


160 


- 


0.99(5) • 10~ 6 


1.06(3) • 10- 6 


165 


- 


1.40(7) • 10- 6 


1.37(3) • 10- 6 


170 


- 


1.81(9) • 10~ 6 


1.53(4) • 10- 6 


175 


- 


2.4(1) • 10~ 6 


1.96(5) • 10- 6 


180 


0.07- 10- 6 


3.1(2) • 10- 6 


2.76(8) • 10- 6 


185 


0.63- 10~ 6 


4.7(2) • 10~ 6 


4.0(2) • 10~ 6 


190 


2.20- 10- 6 


7.4(4) • 10- 6 


6.0(3) • 10- 6 


195 


5.37- 10- 6 


11.3(6) • 10- 6 


9.5(5) • 10~ 6 


200 


10.7- 10~ 6 


17.7(9) • 10- 6 


18(2) • 10" 6 



Table 2: Values of the t — > WZb branching ratio, as a function of the top quark mass, 
obtained in the limit of vanishing Z and W widths (stable particles) and with the convolution 
method. In the last column, the estimate obtained with the decay-chain method from the 
t — > bjiu^uv decay is also shown for comparison. 

3 The t —> WZ b decay rate including the finite- width 
effects 

In this section we discuss the calculation of the t — > WZb decay rate by taking into 
account the large finite-width effects of the W and Z bosons. We present in detail the 
two approaches of the convolution and decay-chain method, and critically compare the 
corresponding estimates for the decay rate. 

3.1 The convolution method 

Since the W and Z bosons are unstable particles with finite widths, their production, in a 
physical process, can be approximately described as the production of real stable particles 
with invariant masses distributed according to a given distribution function. The central 
value and the width of such a distribution are controlled by the physical mass and width 
of the unstable particle. From this point of view, for instance, the t — > WZb decay, if 
occurs for values of the top quark mass smaller than the kinematical threshold, proceeds 
through the production of W and Z bosons with invariant masses close, but smaller than, 
their physical values. 

These observations suggest a very convenient way to take into account the finite-width 
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effects in the t — > WZb decay and, in general, in any process with unstable particles 
produced in the final state. One computes the decay rate of the t — > WZb channel as a 
function of generic values of the W and Z invariant masses. Let us denote this quantity as 
T(k^y, fcf ), where k\y and k\ are the virtualities of the W and Z bosons respectively. The 
t — > WZb decay rate is then obtained by performing a convolution of T(k^, k 2 z ) with the 
two invariant-mass distribution functions 



2 „( m ._ m ,_. /T2 \2 



r(mt—mb) r(mt—rrn, — \/K 7 ) 

Y{t ^WZb)= dk% dk 2 w T(k 2 w , k%) x 

Jo Jo 

x p(tf z ,M z ,r z ) Pi k 2 w ,M w ,T w ) (1) 
The distribution p(fc 2 , M, T) is related to the imaginary part of the gauge boson propagator: 

P (k 2 , M, r) = -hm (- * = 1 ^ (2) 

where M is the mass of the particle and M T(k 2 ) the imaginary part of the vacuum polari- 
zation. In the calculation of T(k 2 ) one may neglect, as a first approximation, the masses of 
the particles flowing in the loops. In this limit, by simple dimensional analysis, one gets: 

MT(k 2 ) = k 2 ^ (3) 

where T = T(M 2 ) is the particle width.0 

In the limit of vanishing width, the distribution function (0) reduces to the delta func- 
tion S(k 2 — M 2 ), and the final particles are constrained on their mass shell. In this limit, 
eqs. (|I|) and (Q) are just a consequence of the optical theorem and the Cutkosky rule to 
extract the imaginary part of the physical amplitude. In general, the procedure outlined 
above, which we refer to as the convolution method, represents the simplest way to take 
into account the finite-width effects in the evaluation of the decay rate, and it has been 
already considered in the literature in different contexts (see for example refs. @, 0]). 

The results obtained for the t —> WZb branching ratio, by using the convolution 
method, are presented in table [| These values are also shown in Fig. together with the 
corresponding values obtained in the stable-particles approximation, i.e. by neglecting the 
finite-width effects. We see that, for a top quark mass around the kinematical threshold of 
176 GeV, the effects of the finite widths increase the branching ratio by orders of magnitude, 

2 Note that the full k 2 dependence of r(/c 2 ) must be taken also into account in the numerator of eq. (||). 
This dependence cancels out the singularity at k 2 — appearing in the propagators and also in the sum 
over the polarization vectors of a particle with invariant mass k 2 : 

In turn, this prescription for the sum over the polarization states is necessary to guarantee the positivity of 
the "off-shell" decay rate T(k 2 v , k 2 z ). We also note that in the calculation of T(k 2 v , k 2 z ), since the invariant 
mass of the W-boson is allowed to vanish in eq. (|l|), we have modified the tree- level expression of the 
W-propagator as indicated in eq. (ph. 
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Figure 2: The t — > WZb branching ratio, as a function of the top quark mass, as obtained 
in the stable particle limit and by taking into account the W and Z finite width effects 
with the convolution method. 



and allow the occurrence of the decay even below the threshold. However, despite the 
large enhancement induced by the widths, for the actual value of the top quark mass, the 
branching ratio is found to be of the order of 10~ 6 , which is probably too small for the 
decay to be visible at the Tevatron Run II or even at the LHC. 

The finite-width effects in the t — > WZ b decay rate have also been studied in ref. . 
By comparing our results with those obtained following the procedure of ref. |2|], we find 
that the values of the decay rate are larger by approximately a factor 3. The reason of 
such a discrepancy will be discussed in detail in the next section. 

3.2 The decay-chain method 

The procedure considered in ref. || to take into account the finite-width effects of the 
gauge bosons in the t — > WZ b decay is based on the study of a process which includes the 
decays of the unstable particles. For illustrative purposes, it is convenient to consider first 
the decay chain: 
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In the calculation of the total rate T(t — > bfiv^vv), the finite- width effects of the inter- 
mediate W and Z bosons are then taken into account by correcting the tree-level expression 
of the gauge boson propagators. In the unitary gauge, by neglecting the fc 2 -dependence of 
the imaginary part of the vacuum polarization, the modified propagator can be written in 
the form g |: 

~9 + 7- ^ /n „ (5) 



k 2 - (m- iT/2) 2 \ * (m-iT/2) 2 / 

where the correction to the longitudinal part of the propagator is necessary to respect 
gauge-invariance. In order to preserve unitarity at high energies, the absorbitive part of 
the triple gauge boson vertex should be also included in the amplitude |10|, [IT| . However, 
in the case of interest, this correction is numerically negligible and for simplicity it has 
been omitted closely following the calculation of ref. || 

The total t — » bfiv^vv rate is then used to compute the ratio: 

r(t -> bfJLUpVU) 



BR(W -> \w^BR{Z -> vv) 

which is assumed to be an estimate of the t — > WZ b decay rate. 

A crucial requirement, for the above procedure to be consistent, is that the full de- 
cay chain, which is considered to estimate T(t — ► WZb), can, at least dominantly, only 
proceed through the initial t — > WZb decay. Otherwise, any direct connection with this 
process is obviously lost. This requirement is only partially satisfied in the case of the 
t — > bfiv^vv decay. Three Feynman diagrams, describing at tree-level the t — > bfiv^vv de- 
cay (with v 7^ Up), are those shown in Fig. [1[ modified to account for the W and Z bosons 
decays into the /zz/ M and vv pairs respectively. These diagrams describe the processes which 
proceed through the initial t — > WZb decay. Besides these graphs, however, other three 
Feynman diagrams also contribute to the decay. They are non-resonant diagrams, similar 
to that of Fig. |3], in which one of the intermediate gauge bosons is produced away from 
its mass-shell. In table § we present the values of the branching ratio obtained by using 
eq. (0). By comparing these values with the t — > WZb branching ratio evaluated by using 
the convolution method we find that the results are in reasonable agreement. It should be 
noticed, however, that this agreement is also a consequence of an accidental cancellation 
between the square of the non-resonant diagrams and the interference among these dia- 
grams and those proceeding through the initial t — > WZ b decay. Ideally, one would like 
to consider in the calculation of the decay rate only these three latter diagrams. However, 
such a prescription does not lead, in general, to a gauge-invariant result. 

The decay-chain method has been considered in ref. to account for the finite-width 
effects of the W and Z bosons in the t — > WZb decay. In that case, the decay rate has 
been evaluated from the ratio: 



BR(W -> iiv ll )BR{Z -> e+e- 



(7) 



which differs from eq. (P) for the choice of the final state. The calculation of T(t — > 
b/iv^e + e~) in eq. (0) requires to take into account a set of 10 Feynman diagrams. The first 
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Figure 3: A Feynman diagram, contributing at tree-level to the t — > bjj,u ll e + e decay, which 
has not been considered in the calculation of ref. 



three diagrams are those proceeding through the initial t — > WZ b decay. This is the only 
contribution which is directly related to T(t — > WZb). Other three diagrams involve the 
radiation of a Z or a W boson from one of the decay product of the first W . These are 
non-resonant diagrams in which the intermediate W is produced away from its mass-shell. 
The last four diagrams are those in which the intermediate Z boson is replaced by a virtual 
photon. In the region of small invariant masses of the final electron pair, these are the 
diagrams which give the dominant contribution to the decay rate. 
Two observations are worth at this point. 

i) In evaluating the Feynman diagrams, the authors of ref. || have missed the contri- 
bution of the non-resonant diagram shown in Fig. [| Although numerically small, in the 
unitary gauge, this contribution is necessary to guarantee a gauge-invariant final result for 
the decay rate. In addition, in the Standard Model, the t — *> b^v^e + e~ decay also receives 
a tree-level contribution from a set of four Feynman diagrams in which the virtual Z is 
replaced by a Higgs boson. Also these diagrams have not been taken into account in ref. . 
We find that, numerically, the contributions of these diagrams is negligible for values of 
the Higgs mass larger than approximately 100 GeV. 

ii) In the theoretical evaluation of the t — > bu.i> fl e + e~ decay rate, all Feynman diagrams 
contributing to the process, at a given order in perturbation theory, must be taken into 
account. However, as discussed before, because of the contribution of the virtual photon, 
with the prescription in eq. (|7|) any direct connection between the t — *> b^u^e + e~ and the 
t — > WZ b decays is lost. To show it clearly, we plot in Fig. [| the differential t — > bjiu^e - 
decay rate, as a function of the invariant mass of the electron pair, as obtained either 
by including in the calculation all Feynman diagrams contributing to the process or only 
the three diagrams which involve the initial t — > WZb decay. Although the latter do 
not form a gauge-invariant set of Feynman diagrams, the comparison is instructive. The 
numerical discrepancy, which is seen between the result of the full calculation and the 
genuine t — > WZ b contribution, it is mainly determined by the contribution of the virtual 
photon in the region of small values of the e + e~ invariant mass. In order to reduce this 
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Figure 4: The differential rate dT / ' d{yr? e J m 2 z ) , for the t — > bfiv^e+e - decay, computed by 
considering all Feynman diagrams contributing to the process or only the three relevant 
diagrams which proceeds through the initial t — > WZ b decay. The top quark mass is fixed 
to m t = l75 GeV. 



contribution, the authors of ref. || have introduced a cut on the values of this mass, 
requiring m ee > 0.8 mz- As can be seen from Fig. ^, this cut is indeed effective to that 
purpose, and in the selected kinematical region the t — > bfj J h , fJj e + e~ decay mainly proceeds 
through the initial t — > WZ b decay. However, with this prescription, the definition itself 
of the t — > WZ b decay rate becomes dependent on the specific choice of the cut. 

As we have shown before, a convenient definition of the t — > WZ b decay rate, which is 
independent of any kinematical cut, is provided by the convolution method. The results 
obtained in this way for the t — > WZb decay rate are roughly 3 times larger than those 
obtained in ref. 0, the main reason being the introduction of the kinematical cut in the 
definition of the decay rate. On the other hand, we emphasize that the calculation of 
ref. is of interest for the evaluation of the signal and the relevant background rates in 
the experimental study of the t —>■ WZ b decay when the b^,v^e + e~ is considered as a final 
state. 

4 Conclusion 

In this paper we have computed the t — > WZ b decay rate, at the leading perturbative order 
in the Standard Model, by including in the calculation the crucial effects of the W and Z 
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widths. We finds that these effects increase the total decay rate by orders of magnitude, 
for a top quark mass of approximately 176 GeV, and allow the decay below the nominal 
threshold. 

In the limit in which the finite-width effects are neglected, the results of the previous 
studies of the t — > WZb decay [@]-|| are in disagreement among each other. We have 
repeated the calculation in this limit, and obtained predictions in agreement with those of 
ref. 0, but with numerical differences with respect to those of ref. |lj and, near threshold, 
of ref. §. 

We have also shown that a practical definition of the t —>■ WZb decay rate can be 
provided even in the presence of large finite-width effects of the W and Z bosons. This 
definition is based on a convolution of the decay rate, computed in the limit of stable 
particles, with two distribution functions for the invariant masses of the gauge bosons. 
The distributions are related to the imaginary part of the gauge boson propagator, and 
are centered around the physical value of the W and Z masses with spread controlled by 
the particle widths. The results obtained in this way for T(t — > WZb) differ, by roughly a 
factor 3, from those obtained following the procedure of ref. @. The main reason of such 
a discrepancy is the introduction of a kinematical cut which is included, in ref. [El, in the 
definition of the decay rate. 

As a final estimate of the t — ► WZ b branching ratio, evaluated by taking into account 
the finite-width effects, we get: 

BR(t -> WZb) ~ 2 x 1(T 6 (8) 

for a top quark mass in the range between 170 and 180 GeV. Eq. (|8|) indicates that, within 
the Standard Model, the branching ratio is quite small so that the observation of the decay 
at the LHC is extremely difficult, if at all possible. Alternatively, the observation of the 
t — > WZ b decay at the LHC with a larger rate would signal the presence of new physics. 
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